function [ sim_struct ] = calculate_synergy_sim( synergy_subject,task,syn_component )

%clear Synergy_Sync;
                    norm = strcat('synergy_c01_',lower(task),'.mat');
                    synergy_norm = load(norm{1,1}); % origine
                    %synergy_norm = load(norm); % test

%% synergy simularity calculation
% muscle vector
[sim_musclevector_matrix] = sim_musclevector_operator(synergy_norm.synergy_subject(syn_component),synergy_subject(syn_component));
[muscle_vector] = muscle_vector_sim(sim_musclevector_matrix);

% time profile 
[time_profile,eigenvalue_timeprofile] =  sim_timeprofile_operator(synergy_norm.synergy_subject(syn_component),synergy_subject(syn_component),muscle_vector);
            
[ sim_struct ] = syn_sim_factor( muscle_vector,time_profile,eigenvalue_timeprofile);

end

